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۰ 0۲081۲6591۷7 70۵-1 متا دتعامصنه21ظ احتمافطمم- رم ۲۲1 اوعا م1 0عامبعاعععه ووه‌باو ماقم 
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1 تاوععءع ۸ 0/21 0 ۳۹۰۳۳۹۹۰۵۵ 1 ۳۱6۷۱۵7۵۵ 1 ۳۹2۵ 





(1 


مت حاحالنت ول عحمصصم[ عوگ موصعم ٩0۵ص‏ ۲۵[ متا اوعع 1 0عتامتع(ععع2 ددع۳٩‏ عصموصمم 1۵7۵1 ۱ مرهج ونطع ص :۳۱۸۵۵96 
وعام‌جصهتهن [16محه عمط عممصتاوع ۲۵ و1 هچره ولا ۵۶ عومجمتای م1 ,060مملاوع۲ص1 و1 قععام‌جصنهعهمر 6216و مه عمرقطه عصمتعوصم-صمه ص۷71 
,16512 نع عط همتصصتامه 4صه حصط‌نم‌علم ۲۷۲ عطا عصلوتا 


ول 1 ,صعلوع) وع۲ ما جصع(جآمترجر لمعصم‌جصیه‌صنط 2 و1 ونطع عصه مصتاوعع ۲0۶ عحصتا عصم1 ۳۵ 2 فعمصتاوومم ص0۴۲ مصتاوع؟ 1۶6 :۱۷۲68000102 
۵5 96 ۲26 )تاو رعتعتالنه؟ ۵۶ عحصتا ع) عمصت0ع۴ ۲۵۶ 0مطعمط 0ع1مم‌مصمورممععع 2 و1 تعظ 1 ,و6وع۲ 116 0عا2ع61عع2 ۲ 50۳۲6۵ وا ممعادا۳0 
2001028 96عط) ۵۲۵۲0۶۵ 1 .عصمتنصمی مصتاجعوه لمحصرمه هماع نی انم ال عمط صعطع منم عفمعععص1 نز وعتصتا غوع عطع ۵۶ 1۵7۵1 
01 م6۵90 [-۲۵6 ]1 0۲0۵۲۵89176 لصا اوعع 1166 260016020640 ووععاو اصهوصومم ۱-16۷۵ 06 ر‌جزهم ونطع صا ,عحصتا مصتصصنا عط ۲۵0۱660 
ممصلصنا عطا ۵۶ وععمصصتاوه ۱۱۱۵000 حصتاحصت‌جهوه عطع مصتصصععع1 هم 60و و1 مصط‌تتممله رس صمتم2لصصل‌م-ممام0ععج] عظ 1[ 
0 860 مایت اصععل1 1۳70 ,جصع(دامجر 0262 عاعام‌جومعصا فطل مصتلك‌صقط ما (مما تهج ۲۵ 2 و1 مصطالءمعل2 فنط 1 .وععاعجصهتن2 
هلاه قصعی عا0112)عحطه‌اصا فطل مصه صمتاعءمرفصا مه ولع1۵7 ووععو 0۶ ع‌حامصربنط رم2ل5 عامزجصنهه عطا 0۶ اعع6 عط فصن مصملم اوعا عظ6 عملمصتایه 

6۰ صمروع عط من 29868860 ع22 1۳001۲0 


عط ععصتا روول .2000 1 عمممصتای ۲۵۷ عظ) عمط 2060فصمجصع و1 16 باع5 01202 ۲62 2 صه بو ممتامابتحصله عط مه 2960 :فمصزل‌صز۳۲ 
و5126 6احرحصصمه المصصه عطع رصمتانل 2 صا .1 میتی ممتاه‌علصصتامون عظع صقمطع 0عصنماداه قه۳ه عععا غصعلمه عمط 2 رل مملنهانت صمتاهملصصتاون 
2 10 1620 ماعمج0ت مصلا0قصعی عتا0112عحصع‌اصاً عععم1 عط عصنه فصمتاععء رفص 0۶ م‌جامصناهر المصصیی معط رولع۲7ع1 ووععاو ۵۶ ع‌حاحصباه للمصمرو عط 

0۲6 611101601 ۰ 


لمع10 مه وا متامجمتنوی طممطعلم بمتمل عصصتهکتا ععلامی من 0عودا و ممتاعءمروصا عنلمنهمر عط) وته‌ورمج ونط م1 تعبیله ۲۷ منوت 
ما عااتوعمم عمط وا ممتاعءم‌عماً متاماصتلصی عط رعاصتحه‌اوممی ۲مماععتاه عمج عصمتم‌تصصنا لمع‌تصطععع و6 ععتل فعمصتاه‌مورمو عباظ ,2۸۵06 
علصاملنه27 ع عمط «محط فانصا اوعع 0۶ وعحصتا عمجم مطا و625 ونطا ما متام روما »مهس معط 6و 0 فقطظ عع‌تاصم‌صصلت‌جنی عظع مه ۳2610166 
6 0 تمطا 4عمصتاوو2 ۲۵ وع‌جرهم ونطا ما رو 56ع۲ عطلا مصتیل مفصلوجر ممصلا صنمعع م2 0ع0عع(آهم م2 فتصتامی عبانم عط «لصم عنم 

مهصا 0۵6۶2 06ممط عوه‌صلاً همعط مه مهو عمرم2طو مه 


اع 1 رصمنءممصنل من رف‌مطعصه لی‌نم‌صآ رصمتنمصصت‌ی 04مطنلمانا م۷( وحصط‌تیم‌عا2 ۳۸۲ :6۳۱۲۵۸08 > 
1 مج صهلوع 


5 7-0۱ (4 1>1160610۵6۴فظ طونم عمطابتظ مصتل هم وعبرمن) 
م 06۲صنا 60 باون مناج ووعععه همه صنه ود 270016 ونط1 .م‌تهم‌وع۱ فصمتامهنن مه مصمتوتنع(۱ 0۶ فتاه[ ,ععقصرعم1 ,1 
1 (۲۳۱۷ ۲ طمتاباحاتاظ وصمصصصصمی 0۳6210176 معط ۵۶ فصمتاتلجم 20 فصصتع] رق)(۵) 
1166۴8665/7/4۰0(۰/ 018.0۲2 حصعومع6 ۲ 6۵1 ۲۲۵://6۲ظ۳) 


تصمیم گیری و نحقیق در عملیات 


دوره ۶ شماره ۴( ۴۰ ۱- 4۵۱۷۰ 





لح مممتوتهع۱۱ که لقصتتاوز 
۵ 6۵۵۲21008 


۲ص حصل- ما0 [ ,۷۳۱۳۱۲۰ 





نوع مقاله: پژوهشی 
برآورد بارامترها با استفاده از الگوریتم ]1۷و بهینه‌سازی طرح آزمون عمر شتاییده تتنش 
ثابت با بارامترهای غبرثایت تحت سانسور فزاینده نوع اول 


9 


گروه ریاضی 9 علوم مهندسی » دانشگاه بین‌المللی امام خمینی (ره) ب مرکز آموزش عالی فتن و مهندسی بوئین‌زهرا بوئین زهرا ایران. 


حکده 


چا 


هدف: در این مقاله آزمون عمر شتابیده تتش ثابت ] مرحله‌ای تحت سانسور فزاینده نوع اول» برای توزیع لومکس در حالتی که هر دو 
ترا ای را مه رس تاره ی هدف از اين مقاله برآورد پارامترهای مدل با استفاده از الگوریتم 1717 و بهینه کردن 


طرح چنین آزمونی است. 


روش‌شناسی پژوهش: اغلب آزمون‌های طول عمر به مدت‌زمان زیادی جهت انجام آزمون نیازمندند و این‌یک مسئله در طراحی آزمون 
است. یک‌راه مواجه با این مشکل, استفاده از آزمون‌های عمر شتابیده می‌باشد. مکانیزم اين آزمون‌ها به این صورت است که با اعمال تنشی 
بیش ازحد طبیعی به محصولات» روند تخریب محصول را افزایش می‌دهند و به‌این‌ترتیب مدت‌زمان انجام آزمون کاهش می‌یابد. آزمون‌های 
عمر شتاییده انواع مختلفی دارند که در این مقاله از مدل تنش ثابت استفاده‌شده است. برای برآورد ماکسیمم درستنمایی پارامترهای محهول 
از الگوریتم /] استفاده‌شده که برای داده‌های مفقودشده روش مناسبی است. همچنین به‌منظور بهینه کردن طرح آزمون موردبررسی از دو 
معیار استفاده شد و تأثیر عواملی همچون حجم نمونه» تعداد سطوح تنش, تعداد بازبینی‌ها و احتمال سانسور میانی بر کارایی آزمون 


موردبررسی قرار گرفت. 


یافته‌ها: بر مبنای مطالعات شبیه‌سازی ویک محموعه داده واقعی نان داده شد که برآورد 1 ارائه‌شده» خوبت ات همجنین نحت معبار 
بهینه‌سازی دوم اه ۱ اول حاصل شد. به‌علاوه برای دست‌یابی به آزمون کاراتر بهتر است از ححم نمونه» تعداد 
سطوح تنش و تعداد بازبینی کمتر و احتمال سانسور میانی بیشتر استفاده کرد. 


اصالت /ارزش‌افزوده علمی: در این بررسی برای جمع‌آوری داده‌های طول عمر از روش بازبینی دوره‌ای استفاده‌شده است. گرچه روش 
بازبینی مداوم یک روش ایده آل است. ولی گاهی اوقات به دلیل محدودیت‌های موحود در آزمون و یا محدودیت بودجه امکان بازبینی 
مداوم وحود ندارد و آزمایشگر ناچار به استفاده از روش بازبینی دوره‌ای است. در این صورت. زمان دقیق شکست محصولات در دسترس 
نیست و فقط تعداد شکست‌ها در زمان‌های خاصی در دسترس می‌باشد. به‌علاوه در این مقاله در نظر گرفته‌شده که هر دو پارامتر شکل و 
انس دزیم توکس رای راک ی با یقت ۰ 


کلبدواژه‌ها: الگوریتم 1۷ تارفتزم دوره‌ای» برآورد ما کسیمم درستنمایی» بهینه‌سازی طرح آزمون روش‌های عددی. 


نو بسنده مسئول 
۴ ۵ ۸۵05۳1۰۳۵1۵۴01 
(4 5 ۱ ۱ 11۱ 


۱- مقدمه 


امروزه با پیشرفت تکنولوژی و قوانین حمایت از مصرف‌کننده. تولیدکنندگان بر این شدند که محصولاتی باکیفیت و قابلیت اطمینان 
بالاتر تولید کنند. این امر. مطالعه طول عمر محصولات باقابلیت اطمینان بالا را تحت شرایط طبیعی استفاده از آن مشکل و گاهی 
غیرممکن می‌سازد. درحالی که اجرای کارکرد محصول در سطح تتش بالاء منجر به کوتاه شدن طول عمر محصول می‌شود. نکته قابل‌تامل 
در مقوله مطالعه طول عمر محصولات این است که. ممکن است شکست یک محصول به مدت‌زمان طولانی مانند جندین سال نیاز 
داشته باشد. یک‌راه مقابله با این مشکل استفاده از آزمون‌های عمر شتابیده! می‌باشد که برای تسریع روند خرابی محصولات و استخراج 
اطلاعات قابلیت اطمینان در شرایط استفاده طبیعی فراهم‌شده است. در محموعه‌ی این آزمون‌ها واحدهای آزمودنی در معرض تنشی 
بالاتر از سطح معمول قرار می‌گیرند تا به این وسیله واحدها سریع‌تر از بین بروند. اطلاعات به‌دست‌آمده از این آزمون‌ها» بر اساس نوع 
رابطه‌ای که بین طول عمر با تتش وحود دارد تجزیه‌وتحلیل می‌شوند. درنهایت نتایج به‌دست‌آمده برای برآورد توزیع طول عمر محصول 
تحت شرایط استفاده طبیعی تعمیم داده می‌شوند. میکرو اسکوبار" (۱۹۹۸) سه فصل کامل را به این مبحث و احرای صحیح چنین 
آزمون‌هایی اختصاص دادند. نلسون" (1۹۹۰) انواع گسترده‌ای از مدل‌های آماری» طرح آزمایش‌ها و روش‌های تحلیلی که در آزمون 
عمر شتابیده مورداستفاده قرار می‌گیرند را ارائه داد. برای حزئیات بیشتر در مورد تحقیقات و مسائل مربوط به آزمون عمر شتابیده به 
بالا کریشنان و هان* (۲۰۰۸) بالا کریشنان و همکاران* (۲۰۰۷) و کتری و کمپس" (۲۰۱۵) مراحعه کنید. 


دو تا از مهم‌ترین انواع آزمون‌های عمر شتاییده؛ تتش گامبه‌گام" و تتش ثابت" هستند. در این مقاله آزمون عمر شتابیده تتش ثابت 
موردبررسی قرار می‌گیرد. در این نوع آزمون» هر نمونه در معرض سطح تنشی ثابت که بالاتر از میزان طبیعی آن می‌باشد» قرار می‌گیرد تا 
زمانی که محصول با شکست مواحه شود و یا با اعمال کردن نوعی سانسور, آزمون به پایان برسد. اگر تنها دو سطح تنش مختلف برای 
آزمون به‌کاربرده شود به آن آزمون عمر شتابیده ساده و اگر از بیش از دو سطح تنش مختلف استفاده شود به آن آزمون عمر شتابیده 


م2 ]۰ ۰ مه ‌ 
جندگانه گفته می‌شود. 


ازآنحاکه داده‌های طول عمر اغلب به‌طور کامل مشاهده نمی‌شوند. حذف واحدهای آزمایش به روشی از قبل برنامه‌ریزی‌شده با به‌صورت 
تصادفی بسیار متداول است. این امر به دلیل زیاد بودن متوسط طول عمر محصولات. کم بودن فرصت اعلام نتایج» محدودیت هزینه و 
عدم قطعیت موجود در نمونه‌های تحت آزمون رخ می‌دهد که به آن سانسور می‌گوييم. یکی از کاراترین روش‌های سانسور» سانسور 
فزاینده نوع "9 است. در این نوع سانسور در پایان هر بازه زمانی از پیش تعیین‌شده» تعداد مشخصی از واحدهای آزمایش سالم» از 
آزمون خارج می‌شوند. دلیل کارایی سانسور فزاینده در این است که واحدهای آزمایشی سالم حذف‌شده از آزمون را می‌توان در 
آزمایش‌های دیگر استفاده کرد. مطالعات زیادی درزمينه استنباط روی آزمون عمر شتابیده تتش ثابت تحت سانسورهای مختلف و برای 


توزیع‌هایی از قبیل نمایی ‏ وایبل " و نمایی تعمیم‌یافته*" انجام‌گرفته است. 


عادع1 110 ۸66۵16۲2160 ۱ 
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براورد بارامترها با استة 


ده از الگوریتم 7و 


سح 


بهینه‌کردن طرح آزمون عمر شتابیده تتش ثاب 


بت 


با پارامترهای غیر ثابت تحت 
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سانسور فزاینده نوع اول 
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۵ در عملیات. دوره 5 شماره ۴۲ ز هت از ۰ ۳۰ صفحه: ۸۵۱۰-۱۱ 


در اغلب قریب به‌اتفاق مطالعات انجام‌گرفته در این زمینه» استنباط بر اساس این فرض انجام‌شده که فقط یکی از پارامترهای مدل با تتش 
اعمال‌شده بر واحدها مرتبط است و سایر پارامترها در تمام مراحل آزمایش ثابت می‌مانند. به‌عنوان‌مثال به پاسکال" (۲۰۰۸ میترا و 
همکاران" (۲۰۱۳). گوان و همکاران " (۲۰۱6) و فن و هسو* (۲۰۱) مراحعه نمایید. بااین‌وحود» در عمل مشاهده می‌شود که این فرض 
ممکن است مناسب نباشد و مدل با یک پارامتر ثابت که وابسته به سطح تتش نیست» نمی تواند برازش مناسبی برای داده‌های واقعی باشد. 
در اغلب مقالات که برای طول عمر واحدهای آزمون از توزیع طول عمری با بیش از یک پارامتر استفاده‌شده است؛ معمولاً فروض شده که 
پارامتر مقیاس * وابسته به سطح تنش است ولی پارامترهای مکان" و شکل " برای تمامی سطوح تدش ثابت هستند. درصورتی که بررسی‌هایی 
وحود دارد که این مطلب را نقض می‌کنند. 


به‌عنوان‌مثال» لی و همکاران" (۱۹۸۹) و هیرحیست و همکاران" (۱۹۸۹) دریافتند که پارامتر شکل توزیع وایبل به ترتیب:فاثر از تشن 
میدان الکتریکی برای دی‌الکتریک‌ها و تتش دما در آزمایش خازن‌ها می‌باشد. به‌علاوه» نلسون (۱۹۸) و بوی‌کو و گرلاج " (۱۹۸۹) 
مشاهده کردند که در تحزیه‌وتحلیل داده‌های واقعی زمانی که توزیع مکان-مقیاس برای طول عمر واحدها استفاده‌شده. هر دو پارامتر مکان 
و مقیاس به تتش بستگی دارند. وانگ" (۲۰۱۷) خانواده توزیع‌های بریده‌شده"" دو پارامتری و وانگ (۸ توزیع وایبل را در حالتی که 
هر دو پارامتر آن‌ها با تغییر سطح تنش غیرثابت می‌باشند برای آزمون عمر شتابیده با دو سطح تنش (ساده) موردبررسی قراردادند. ازآنحاکه 
آزمون شتابیده ساده» به لحاظ محاسباتی ساده‌تر است ولی آزمون شتابیده چندگانه کارایی بیشتری دارد» در این مقاله قصد داریم آزمون عمر 
شتابیده تتش ثابت را در حالت چندگانه برای توزیع لومکس " و تحت سانسور فزاینده نوع اول مورد تجزیه‌وتحلیل قرار دهیم. همچنین در 
نظر می‌گیریم هر دو پارامتر این توزیع متأثر از سطح تدش هستند. توزیع لومکس برای اولین بار توسط لومکس ( ۱۹۵) معرفی شد. تعداد 
کمی از محققان این توزیع را برای آزمون‌های عمر شتابیده استفاده کردند. به‌عنوان‌مثال حسن و همکاران* (۲۰۱7) و چاندرا و خان*" 
(۲۰۱۸) با در نظر گرفتن توزیع لومکس به بهینه کردن آزمون عمر شتابیده تتش گام‌به‌گام تحت معیارها وفرضیات مختلف پرداختند. 


به‌طورکلی دو حالت بازرسی در آزمون‌های عمر شتابیده به‌کاربرده می‌شود. یکی بازرسی مستمر" " و دیگری بازرسی دوره‌ای"". بازرسی 
مستمر معمول‌ترین روش بازرسی است که در آن زمان دقیق خرابی موجود است؛ اما در بازرسی دوره‌ای» زمان دقیق خرابی در دسترس 
تست دس سیر کح هلکسا کی تم مس هناسفا 1 
روش استفاده کرد. در این صورت به‌حای زمان دقیق شکست محصولات. تعداد شکست محصولات در بازه‌های زمانی از پیش تعیین‌شده 
ثبت می‌شود. از مزایای استفاده از بازبینی دوره‌ای می‌توان به‌سادگی اعمال این روش ارزان و در دسترس بودن آن نام برد. در این مقاله 
روش بازرسی دوره‌ای را برای جمع‌آوری اطلاعات تحت آزمون تتش ابت مورداستفاده قرار می‌دهیم. برآورد پارامترهای مدل با استفاده از 
الگوریتم اتتظار تیه که به‌اختصار با نماد 71 نمایش داده می‌شود. ارائه می‌شود. یکی از مهم‌ترین اقدامات در مواجه با آزمون‌های 
طول عمر بهینه کردن طرح این آزمون‌ها است. در اینجا هدف از بهینه کردن طرح آزمون انتخاب بهترین زمان برای خاتمه آزمون است. 
البته این زمان بهینه» متأثر از انتخاب معیار مناسب حهت بهینه‌سازی است. در این مقاله از دو معیار برای بهینه‌سازی استفاده می‌شود که 
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هر دو بر یایه ماتریس اطلاع فیشر می‌باشند. درنهایت» تأثیر احتمال سانسور میانی» ححم نمونه» تعداد سطوح تنش و شاخص‌های دیگر 
در طراحی آزمون بهینه موردمطالعه قرار می‌گیرد. 


در بخش ۲ مقاله حاضر به توصیف مدل موردبررسی و فرمول‌بندی آن پرداخته می‌شود. در بخش ۲ الگوریتم 7711 برای به دست آوردن 
برآورد ماکسیمم درستنمایی پارامترهای مجهول ارائه می‌شود. بهینه‌سازی طرح آزمون موردبررسی با بهره‌گیری از دو معیار مختلف در 
بخش ؟ موردبررسی قرار می‌گیرد که به دلیل این که هر دو معیار بر مبنای ماتریس اطلاع فیشر هستند» این ماتریس نیز در همین بخش 
ارائه می‌گردد. بخش ۵ نتایج مطالعات عددی انجام‌گرفته» در راستای مطالب تئوری بیان‌شده و همچنین بررسی عوامل موثر بر طراحی 
آزمون بهینه را ارائه می‌دهد. بخش ۶ به بررسی یک محموعه داده واقعی می‌پردازد. نهایتاً نتایج به‌دست‌آمده از این بررسی در بخش ۷ 


بیان می‌شود. 

۲- فرضیات اساسی و روش انجام آزمون 
۲-۱ - فرضیات اساسی 

فرضیات زیر در این مقاله در نظر گرفته می‌شوند: 


طول عمر واحدهای آزمون تحت سطح تنش نرمال مک و سطوح تنش شتابیده ٩,‏ > ...> رک > رگد از توزیع لومکس با تابع توزیع زیر 


تبعیت می‌کند: 
)۱( و وش ول و0 < 1 ,ما < 1 , + - 1 < (۸) :7 


برای بیان رابطه بین تتش و طول عمر محصولات از مدل لگ خطی برای هر دو پارامتر شکل و مقیاس به‌صورت زیر استفاده می‌شود: 


)۲( .رگرب + وه < ب2 108 
۳( ۰ رصرتر + ور < 0 102 


به ازای 6و 2 ,0 < 7 که 0۷ ۰ ,0۶ ) م0 و ر صرایبت محهولی هستند که باید برآورد شوند. 
۲-۲- روش انجام آزمون 


۱ آزمون عمر شتابیده تتش ثابت» در / سطح مختلف برای نمونه‌ای به ححم از محصولات به‌طور هم‌زمان انجام می‌شود. 
3 

به این صورت که در سطح تنش شتابیده رگ که /,...,0,,2< ۰ تعداد ,70 واحد از محصولات تحت آزمون قرار می‌گیرند. 

51 در سطح تنش ]ام تعداد 71 بازرسی دوره‌ای در زمان‌های از پیش تعیین‌شده روآ و--و زگ و رز نورب فین گرد که نآ زمان خاتمه آزمون در 
سطح تیش ام می‌باشد. 

۳ _در سطح تتش ,گ» تعداد واحدهایی که در بازه رز ,ر-ز) با شکست مواجه می‌شوند. ثبت می‌شوند که آن را با نماد ررتد» نمایش می‌دهیم که 
0 1 و وه صولولا ای ور اک از : 

4 سانسور فزاینده نوع اول را برای صرفه‌جویی در زمان و هزینه روی واحدهای آزمون اعمال می‌کنيم. به این صورت که در زمان ررل» ,ره تا از 
واحدهای سالم به‌طور تصادفی از آزمون حذف می‌شوند؛ به‌طور مشابه در زمان مر رت تا از واحدهای سالم به‌طور تصادفی کنار گذاشته 


می‌شوند و این روند به همین ترتیب ادامه می‌یابد تا زمانی که به آخرین زمان بازرسی که معادل بازمان خاتمه آزمون در سطح تنش رگ می‌باشد 
[- :171 171 


سم 1 
یعنی نآ می‌رسیم؛ 5۳ ات زمان تمامی واحدهای سالم باقی مانده از ازمون حدف می‌شوند؛ یعنی رز 2 زد مد س_ :11 -- رز ۲ 
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۵ در عملیات. دوره س شماره ۴۲ ز هت از ۰ ۳۰ صفحه: ۸۵۱۰-۱۱ 


۳- برآورد پارامترهای مجهول مدل 


روش‌های مختلفی برای تخمین داده‌های سانسور شده وحود دارد. در اینحا از روش برآورد ماکسیمم کشت تما نی به دلیل سادگی 
محاسبات و برخورداری از توزیع مجانبی نرمال و تخمین واریانس مجانبی برآوردگرها؛ استفاده می‌شود (نلسون, ۱۹۹۰). تابع درستنمای 
تحت سانسور فزاینده نوع اول با حذف‌های تصادفی؛ برای مشاهدات زد #,..,1-0,1,2 و ,...,1,2- زر و نسبت به پارامترهای 
وب ی )0 باصورت زیر ماهس رده 


1,60 (1)6۱6(۳)6(۰ (۴( 


می‌توان نوشت: 
۳۵ 2 6 ارت (وهه رنه اجه)ط( نا جه)ط(ره) - ( زه) ۲ 
۱ ۱ كِ 7 
(۵) / ۲ 
یم رش اج 1( ۳ ,2 5 رم " ۱ [ 1 س_ 
ز [< [ [-1 


][- 170 17 
به یاد داریم که 6 3 وود و2 رورس اه 
1-1 [-] ۱ 
از طرف دیگ با توحه به اینکه مشاهدات ند ]» از توزیم چندجمله‌ای با پارامترهای ,۵ و رم پیروی می‌کنند که بر احتمال شکست 


محصول در بازه زمانی (ز ۰ -ز) می‌باشد. پس 


۳ 
ِ ِ 7 / 1 ۰۰ ۰۰ 

(۶ . "7 "رم ] ]1 ره !زد ۱11۷:۱111 -(1)0۱۶ 
ِ [-ز 1 آ-ز [-| [-| 


که در رابطه (۶ داریم: 





ی پل 
)۷( ۱ | + ]دراو وه 





7 2 .1 
۸ ۳ -۳)۶,۱ -7- 
۱ + ( )1۳ 1 یط 
با جایگذاری روابط (۷) و (۸) در رابطه (۶) و نهایتاً با حایگذاری روابط (۵) و (۶) در رابطه (۴) تابع درستنمایی به‌صورت زیر 
محاسبه می‌شود: 
نزن ۹ 2 .1 4 ریگ ز1 _6/ 
۹( 1 رب + | |[ | 1/0۶ 
0 0 0 


جمتاحصصتا و۳ 0ممطتای121 صیاجورن‌بع]۱۷ . !1 


اکنون با استفاده از روابط -خطی که د بین پارامترهای نوزیع لومکس و سطح تنش به فرم روابط ۲۱( و (۳) در نظر گرفته‌شده بود» 
می‌توان لگاریتم تابع درستنمایی ارائه‌شده در رابطه )٩(‏ را نسبت به پارامترهای (,/.ر00.0,۰/5)< ۶ به‌صورت زیر نوشت: 


( 6۱۵۵ | روروز) ما-(وهد رن | وم)گ( رها مه) ۴( )۳ - ( زه) ۲ 





۱۰( ۳ بٌ_ فا ] 
ك ِ ۳( - 7) ۳ ۸ ِ ا]] ِ ۱ 11 - 
زز6 [<-[ 1 5۷۵ 


واضح است که با به دست آوردن برآورد پارامترهای 7 با توحه به روابط (۲) و (۳) می‌توان برآورد پارامترهای را نیز محاسبه کرد. 
برای به دست آوردن برآورد ماکسیمم درستنمایی پارامترهاء باید از رابطه (۱۰) نسبت به پارامترهای م۰6 ,» م5 و ,۸ مشتق بگیریم 


و آن‌ها را برابر صفر قرار دهیم: 


۳ 

2 

(ع این 0 

| وه + درک ۲0 19 00*0٩‏ زه )9 - 2 

(۱۱ ار سر 2 
نمی رورم "سم و 1 

7 << [<-[ 9 

زهر6 + 6+٩‏ ۲ این 0 

| و + 9 » رک رنه زد ۷ ۱ 

[<[ <ا نه 

۱۲( 
۱ ۳ ۳ )وه ۹ درد( - ( ,کر + رق) و رک بثِ ۱ 

ً ۳ 5" 

مرن ود - (>د,() 0 

9 نف كت [-ز | متا0 3 

۱۳( ۱ ۱ 
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همان طور که در روابط (۱۱-(۱۴) مشخص است. به دست آوردن برآورد ماکسیمم درستنمایی پارامترهای م۵ ربه» رت و ,8 به‌صورت 
تحلیلی امکان‌پذیر نیست؛ ازاین‌رو باید از روش‌های عددی استفاده کنیم. به‌عنوان‌مثال می‌توان به مقالات خواحه‌زاده و همکاران" (۲۰۲۰) 
برنیانی‌فر و همکاران" (۲۰۱۸) و خوش‌فطرت و حسین زاده لطفی " (۲۰۱۴) اشاره کرد. در اغلب مقالات» روش نیوتن-رافسون * برای به دست 
آوردن بر آورد ماکسیمم درستنمایی در چنین شرایطی مورداستفاده قرارگرفته است. به‌عنوان‌مثال می‌توان به مقاله برجمی و همکاران" (۲۰۱۹) اشاره 


بت تحت 
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۵ در عملیات. دوره 5 شماره ۴۲ هت از ۰ ۳۰ صفحه: ۸۵۱۰-۱۱ 


کرد. برای به کاربردن روش نیوتن-رافسون» مشتقات مرتبه دوم تابع لگاریتم درستنمایی در هر تکرار موردنیاز است که اين امر ممکن است 
گاهی اوقات پیچیده باشد؛ بنابراین برای حلوگیری از این مشکل» از الگوریتم 171 در این مقاله استفاده می‌کنيم. همچنین لیتل و روبین" 
(۲۰۰۷) اشاره کردند که گرچه الگوریتم /77 سرعت همگرایی کمی دارد اما در مقایسه با روش نیوتن-رافسون قابل اعتمادتر است خصوصاً 
زمانی که تعداد زیادی از داده‌ها نافقص باشند. به الگوریتم 7 در بخش بعدی اشاره می‌شود. 


۳-۱- الکوریتم 11 


الگوریتم 711 یک روش معمول ترا ماکسیمم کردن توابع درستنمایی بیچیده در برخورد با داده‌های ناقص اتف ا: فرای اسان پیشتر با این 
الگوریتم می‌توانید به رشیدی نژاد و نواب پور" (۲۰۱۰) مراحعه کنید. برای استفاده از الگوریتم 7 ما دودسته داده داریم: دسته اول. 
محموعه داده‌های مشاهده شده (داده‌های ناقص) و دسته دوم محموعه داده‌های غیرقابل مشاهده (داده‌های کامل). در واقع محموعه 
داده‌های ناقص مها مشاهده می‌شوند و محموعه داده‌های کامل همراه با تعدادی داده گم‌شده می‌باشد. ما قزر این مقاله زمان سس 
دقیق مولفه‌ها ر به‌عنوان داده‌های گم‌شده» در نظر می‌گيريم. به‌عبارت‌دیگر در نظر می‌گیریم که مولفه‌ها نحت بازرسی مداوم و بدون اعمال 


حقیقت دنباله مشاهدات |ز] که قبلاً معرفی کرده بودیم. تابعی قابل مشاهده از داده‌های کامل است. لفظ دادهای ناقص در مورد داده‌های 
گم‌شده به کار می‌رود اما به هر داده‌ای که هرگز قابل‌مشاهده نباشد نیز اطلاق می‌شود. 


در الگوریتم /۷/ روش به این صورت ینت ) کف اتنت! پاید ثرا داده‌های کامل که پر کی از داده‌های مشاهده‌شده و داده‌های گم‌شده می‌باشند» 
لگاریتم تابع درستنمایی را محاسبه کنیم که آن را با نماد (2),رر./ نمایش می‌دهيم. الگوریتم ۸ با استفاده از لگاریتم تابع درستنمایی 


۱ ۲ ۱ ۲ شاک( 
داده‌های کامل طی فرایندی تکراری به طور غیرمستقیم» فاد لد درست: ی و( اج ۹ 


را حل می‌کند. فرض کنید (7 مقداری 
آغازین برای 7 باشد؛ در تکرار اول گام 7 عبارت است از محاسبه (ت«ل(),رمی/)2 و گام /۸ شامل محاسبه ۱ از طریق ماکسیمم کردن 
(د()برم/) 17 . گام‌های 17 و ۸۸ مجدد تکرار می‌شوند. البته این بار به‌جای ( از !7 استفاده می‌شود. در تکرار (+:) ام» گام‌های 


5 )ریم )1 باید محاسبه شود. 


5 انتحاب مقدار ۳ به‌طور یکه تابع )1 را مأکسیم مکند. 
گام‌های ۶ و ۸۸ تا زمانی تکرار می‌شوند که تفاضل مقدار کوجک دلخواهی را که نشان‌دهنده همگرایی دنباله است» اختیار کند. 
از مشاهدات ثبت‌شده و مشاهدات سانسور شده می‌باشد؛ بنایراین ما این گونه در نظر می‌گیریم که بیان‌گر طول عمر واحدهایی هستند که 


در بازه زمانی با 0 مواحه شده‌اند و بیانگر طول عمر واحدهایی تفت که در زمان بازرسی 7 به‌تصادف از آزمون کاز 9 میی 


شوند؛ بنابراین تحت سانسور فزاینده نوع اول نمونه داده‌های کامل را می‌توان به‌صورت نوشت که. 


درنتبحه با توحه به نکات ذکرشده و روابط (۲) و (۳) لگاریتم تابع درستنمایی داده‌های کامل به‌صورت زير محاسبه می‌شود: 


/ / 
(رکر + وک ) ۷ - (رفره + مه) ۷ ( ۳ ,)ریک 
1 ] 


(۱۵( 
زفرا- ملسم * ك زهرثا- متا | 1 00 
۱ [ ۳ +1)عها ‏ +[ + ِِِپِ )۵ 2 7+ ۳ ۰( - 
[-/ -ر |1<[ [< 


صزجای ۲۲ 0صیه 110016 ۱ 


۱۵۷۷۵۵۵0۵۴ 0صح 20[مصلل تطوم ۶ 


برای احرای مرحله ل] داریم: 


ّ 3 3 
بر + + ر + 7 یز 
5-2 وب (7+ )وا کیب ۵+ )۷ --(۱ اد |[ ان 





4 4 [-<1 
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و 
4 
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اکنون برای اجرای مرحله ۸۸ به مشتقات جزنی رابطه (۱۶) نیاز داریم: 


1 [+۱] 
ریب ( ۳و ِ / ۱ از | بر ۱ 0 هگ 9۳ 


0 (#ر < 


۱ 5 
)ره )رمع 
۳۳ 
کي »5 5 و ۳ تسا اسف سب 
۱۸ 7 ۳ ۲ 
6 ( ام 0 7 9 "1 ۳ ۳۳ ,( ( یم رت )۴ 2 
۳2 
0 ۰ زد ۱ 7 7 ی 
(۱۰) 9 ۲ 77۳ 


2 
0) 


رهز 


که در روابط (۱۷)-(۲۰) داریم: 


۳ 1 [< وتا 
[ 8 0۳ ِ ش () 10 ۲ ره) ۲ +۳۱ ۳ () ]وت را 


رز > ره |( + ۷8 )۱۵ ]2 < (,۵۰/(۰.مه) رت 


۷ | ند | بز ,رز و 
)" ر +( ون ,+ مکی زک ۳۹ 


۲۱0( ۱(۶)2(۵۴+ ۳ ره )ها ۱ 7 


| < || 1+ نطر وی )10۵ سس ( رگا متا ,ره منت ) #9 


(۳۲ ۲ )+ نمی ]یک ۳ 


زه | 9 +1 
0 


درز ] زا ۱ 1+ 0۵ رح - ( رک ,مک +ره ,مت ) و 








(۲۳( ِ زا 1 
11 ۳۳ 
9 ندر 0 ی 7 أ 9 7 ۱ 
۳ ات ا ۳ 
0 0 
ً ۳ 
مسب و رواب - / / و من 2 
7 7 / [ + 0/۸۹ رج ۱ 0۰۱۳۰۵ + 
(۲۴( و9 
۷ 11 / 
7 7 ِ 
۱ ۱ نکر-0 بح 71 3 
3 | +1 


/ 
4 


۰ 
۰ 
فب 


که (2) ۶ تابع احتمال توزیع لومکس می‌باشد؛ به‌عبارت‌دیگر مشتق رابطه (۱) است. 


۴- بهینه کردن طرح آزمون 


حکمی‌پور / ت 


همان‌طور که قبلاً اشاره شد یکی از اهداف این مقاله بهینه کردن طرح آزمون تتش ثابت است. منظور از بهینه‌سازی طرح آزمون تدش ابت» 
تعیین مدت‌زمان بهینه اعمال هر سطح تنش است. معیارهای مختلفی برای بهینه‌سازی مورداستفاده قرارمی‌گیرند که همه آن‌ها مبتنی بر 
ماتریس اطلاع فیشر می‌باشند. در اين مقاله از دو معیار برای بهینه‌سازی طرح آزمون استفاده می‌شود. معیار اول عبارت است از به حداقل 
رساندن واریانس مجانبی برآوردگر میانگین طول عمر محصول و معیار دوم نیز شامل به حداقل رساندن واریانس مجانبی برآوردگر چندک 
مرتبه 0 ام می‌باشد. 


گیری و ند 


ف 
مه 


۳-۱ - معیار بهینه‌سازی اول 


این معیار از متداول‌ترین معبارهای بهینه‌سازی در مطالعات طول عمر اشتت) و زمانی مورداستفاده قرار قم کر که میانگین طول عمر 
محصولات اهمیت وا 3 باشد. در این بحش قصد داریم واریانس محانبی برآوردگر میانگین طول عمر محصولات در شرایط استفاده 


معمول را به حداقل برسانیم. 


میانگین طول عمر محصولاتی که از توزیع لومکنن عبت می کننك» در شرایط استفاده معمول به‌صورت زیر می‌باشد: 





۵ در عملیات. دوره س شماره ۴۲ هت از ۰ ۳۰ صفحه: ۸۵۱۰-۱۱ 


وتا + وت 
[- ۲,50مللی ‏ 1- م2 
برای سادگی محاسبات از لگاریتم برآوردگر میانگین طول عمر محصول استفاده می‌کنیم: 
(۲۶) |1- ۷و ) عه( - ور + م6 > ( و ) 102 


واریانس محانبی م/ 102 را با نماد ( رت 102) ۸۲۰۲ نمایش می‌دهیم و برای محاسبه آن از روش دلتا استفاده می‌کنیم: 
۳۲۷( ۸( ۸۷۵۲)۱02 - 2 


که ۶ عبارت است از برآورد ماتریس اطلاع فیشر و رآ بردار مشتقات مرتبه اول از ,108 نسبت به پارامترهای م۵ ر۵» مت ور 


می‌باشد. در عمل مقادیر 0,۵ 4) - 2 بر مبنای یک تحربه قبلی و یا بر اساس نتیجه آزمون مقدماتی به دست می‌آیند. طبق 
توضیحات ارائه‌شده ,77 به‌صورت زیر محاسبه می‌شود: 


00-0 00-0 
0 نم 0 نم 


۳۸( ماگ7 
0 7 00 - 00 1 


اکنون به محاسبه ماتریس اطلاع فیشر! می‌پردازيم . همان‌طور که می‌دانیم» از کاربردهای ماتریس اطلاع فیشر تعیین وازنانس بر آوزدگر 
و رفتار محانبی برآوردگر ماکسیمم درستنمایی است. ماتریس اطلاع فیشر» دارای نقش کلیدی در برآورد پارامترهای مدل می‌باشد و 
توصیف‌کننده میزان اطلاع داده‌ها درباره پارامترهای مجهول است. به‌علاوه معکوس این ماتریس برآوردی برای ماتریس واریانس- 
کواریانس مجانبی برآوردگرهای پارامترها می‌باشد. برای به دست آوردن ماتریس اطلاع فیشر باید از منفی مشتقات جزئی مرتبه دوم 
لگاریتم تابع درستنمایی نسبت به پارامترهاء امید ریاضی گرفت. در اینجا ما ماتریس اطلاع فیشر را با نماد * نمایش می‌دهيم. با اين 


توضیحات. درایه‌های ماتریس اطلاع فیشر به‌صورت زير محاسبه می‌شوند: 


111 1 
| + نرق +20 19 زظ 0+6 >( زه - رد)2 < <- ر رت 
(۲۹) ٍِِ« 
1 6 17 ] 
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[ز | ِ 


۷2( ممتاجصمکو مور[ 





۰۷۹ 


سم 


براورد پارامترها با استفا 


ده از الگوریتم 12۷۲ و بهینه‌کردن طر 


سح 


ح آزمون عمر شتابیده تتش ثابت ؛ 


با با 


ک‌ 


رامترهای غیر ثا؛ 


بت تحت 


4 


سانسور فزاینده نوع اول 





(۳۴( : ۳ 
+00 ۵ 0+05 ۹ 
كٍِِ ک(زه)تر ور ۲ ک(رداظ + 
هرت +۲0 
[< [ << [ 1 كُ [< ] 
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3 0۵40 2 
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۸۰ :03-0 2 ۱ ح 3 0/308 2 ۱ 
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۱ 0۵ بو 
| + ۱ 0 ۴ [-ز -ز 
(۳۶( 
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| 2 ۱ زرف +60 1 
2 بو + تفت ۳1-7 / 6 [-ر از 
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: زد رت 2-2 زرفرت +20 1 
0 0 ۱ 
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ِ + تفا زا - ۹( [-ز | 
۱ (۲۳۷) 
ِ زکرظ 2+ وتاکی در مکی 0 ] 
42 شش مر مس ار "۳ رک( رد2 << + 
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زگرک 2+ وتر2 زفرگر +50 1 1 
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۱ ۸۵ 0 ِ_ 
(۳۸) 


زگرک 0+2 ۳ ره + نی 


1 719 
+ 1) تک‎ 0*٩ 


۱2 میم ار مر مر 
وه + تفت ز- ۳1 یت [-ز -ز 


که در روابط (۲۹)-(۳۸» دا و (زهات به ازای #,.,0,1,2< 1 و ...1,2 طبق مقاله بودهیراحا و پرادهان! (۲۰۲۰) از 


تا دوره س شماره ۲ ز هت از ۰-۰۰( صفحه: ۰-۲۱ ۷ ۵ 


روابط (۳۹) و (4۰) به دست می‌آیند. 
[-ز 
(۳۹( (,0-72] [ 0 10 [ رد)1 
[<] 
[-[ 
(۳۰( ( 7 )| [زنه ,2 | )17 
[<] 


که‌رر] و۲ در روابط )۷۲( و )۸( ارائه‌شده‌اند. 


مم0 ۳۵ 0صه وزهتتط0 با ! 


۴-۲- معیار بهینه‌سازی دوم 


معیار دیگری که در اینحا برای به دست آوردن طرح آزمون بهینه از آن استفاده‌شده. مینیمم کردن واریانس محانبی برآوردگر ماکسیمم 
درستنمایی طول عمر صدکی تحت شرایط استفاده معمول است. به این منظور ابتدا چندک مرتبه 4 ام طول عمر محصول. تحت شرایط 
عملیاتی معمول» م5 » ارائه می‌گردد که آن را با نماد نمایش می‌دهیم: 


00-0 سم ۳ 
۳۱( - ۱۱-۸ رای - و 


برای سادگی محاسبات از لگاریتم چندک مرتبه 4 ام طول عمر. استفاده می‌کنیم وبا استفاده از روش دلتا واریانس محانبی برآوردگر 
ماکسیمم درستنمایی و102 را به‌صورت زیر محاسبه می‌کنیم: 


(۴۲) .۶ ول - ( راعه1] :۵۷۵ - و 


که ۶ برآورد ماتریس اطلاع فیشر و ۰ بردار مشتقات مرتبه اول از 1081 نسبت به پارامترهای م9 9 0 و / می‌باشد؛ بنابراین 


ماتریس 77 به‌صورت زیر به دست می‌آید: 


00-0 


(108)1-6 (و -1) 2 9 
: 5 ۳ 2 
۳ 
ِ« 0ب درش-م۵- 
۳ ( -1) 108 ( و -1) 10 0 59 
ٍٍِِِِ 00-0 -م.- 
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۵- مطالعات شبیه‌سازی‌شده 


هدف اصلی از این مطالعه بررسی عددی برای نشان دادن نتایج نظری بحث شده هست؛ که این بررسی هم شامل مسئله برآورد 
پارامترهای مدل با استفاده از روش الگوریتم 7711 و هم شامل طراحی آزمون بهینه می‌شود. 


برای این‌که بهتر بتوانیم روند طرح آزمون بهینه را به نمایش بگذاریم از مدل استاندارد استفاده می‌کنیم. تلو ار استت داوزم.حلنت 
مسئله به‌حای استفاده از سطوح تنش اصلی. می‌توان از سطوح تنش استانداردشده که به‌صورت زیر محاسیه می‌شوند استفاده کرد: 


0 - ری 5 


(۴۴ ۱7 و ت۱3 
رک - رگ 


که در این صورت 0 5 و ار روز 


برای سادگی محاسبات می‌توان تعداد بازرسی در سطوح تنش مختلف را با یکدیگر یکسان در نظر گرفت؛ در این صورت داریم: 
0 -.- و191 - ,۰1 به‌علاوه می‌توان در نظر گرفت زمان‌های بازرسی در تمامی سطوح تنش با یکدیگر برابر هستند؛ به‌عبارت‌دیگر 
پر ...۰ رواع ررا به آزای 2,..,۷,/<. در این صورت زمان خاتمه تمامی سطح تنتش با یکدیگر برابر می‌باشند. یعنی 


27 1 < .2۰.۰ روا بر . به‌عبارت‌دیگر 7 زمان خاتمه آزمون است. 


برای ارزیابی برآوردگر ماکسیمم درستنمایی که توسط الگوریتم 21 به دست آمد و در بخش ۳ موردبحث قرار گرفت» از سه ترکیب داده 
مختلف استفاده می‌شود. مثال ۱ که آزمون عمر شتابیده تنش ثابت در دو سطح مختلف اجرا می‌شود و مقادیر اولیه‌ای آن عبارت‌اند از: 





۸۱ 


سم 


براورد پارامترها با استفا 


ده از الگوریتم 12۷۲ و بهینه‌کردن طر 


سح 


رح آزمون عمر شتابید 


۵ ند تدش ثابت 


با پا 


رامترهای غیر ثا؛ 


بت تحت 
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سانسور فزاینده نوع اول 
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حکمی‌پور / ت 


یم‌گیری و تحقیق در 


مب 


تا دوره س شماره ۲ ز هت از ۰-۰۰( صفحه: ۰-۲۱ ۷ ۵ 


2< 1 02 رگ [ع وگ 19-6 40 عرص 50ع ره 03-< مه 05< 9 ۰-04 0.6 ,0 و زمان‌های بازرسی را ۱, ۰۰ 
۰۲ ۰۳ ۰۰۰۴ ۰/۶ در نظر می‌گيريم. مثال ۲» شامل آزمون عمر شتابیده تتش ثابت در سه سطح مختلف است و مقادیر اولیه 
ای آن عبارت‌اند از: 3< 1 [0< رگ 03ع رگ [ع وگ 924 30ع ری 40ع ره دعب 01ع6» ۰0,۶02 8-02 
۰ 2-3 ,8 و زمان‌های بازرسی را ۰۰۱ ۰۰,۲ ۰/۳ و ۰/۴ در نظر می‌گيريم. مثال ۳ نیز شامل آزمون عمر شتابیده تتش ثابت در چهار 
سطح مختلف است و مقادیر اولیه‌ای آن عبارت‌اند از: 4< ۰ 5-202 ۰ 204 وق 206 وک [ع پگ ۰7-6 30 ره 240 را 
» 50< و 60پ 6-12 6[ ره » 7.1 27.6 ,8 و زمان‌های بازرسی را ۸۱ ۰۰ ۰۰,۲ ۰,۳ ۵۰۰۴ ۰و ۱/۶ در نظر 


می‌گیریم. 


همچنین در هر بازبینی در نظر می‌گیریم تعداد واحدهای سانسور شده از توزیع دوحمله‌ای با احتمالات ثابت 6۰ 0۰ و 2۲۰ بیروی 
وب کال واضح است که زمانی که احتمال سانسور میانی معادل ۵۰ است؛ درواقع ما با سانسور نوع اول سروکار داریم. 


تحت این تنظیمات و با ده هزار بار تکرار الگوریتم» نتایج مطالعات شبیه‌سازی در جدول‌های ۱ تا ۳ به ازای احتمال سانسورهای مختلف» 
ارائه‌شده است. اریبی! هر برآوردگر» میانگین توان دوم خطا" برای هر یک از بارامترها که به‌اختصار با نماد ۸/5 نمایش داده می‌شود و 
خطای نسبی " آن‌ها که به‌اختصار با 02 نمایش داده می‌شود. محاسبه‌شده و ارائه‌شده است. همان‌طور که در حدول‌های ۱ تا ۲ مشاهده 
می‌کنیم. مقدار 1/57 به ازای تمامی برآوردگرها مقدار بسیار کوجکی است و این نشان از خوب بودن بر آوردگر ارائه‌شده دارد؛ از طرف دیگر 


واقعی پارامتر مجهول فاصله می‌گيريم. 


جدول ۱ - برآوردگر ,]1۷ بارامترها به همراه اریبی ۷15۳ و 1317 برای داده‌های مثال ۱. 
۰ 62۵100016 ]0 0202 ۲۵۲ ,۱۳ 200 به۷۲ تم ۱۳۱6 ماه همه هن عمط ۵۶ ]۱۷ 1۰ 12016 


71 ,1۱۷1۳ 35( 9۳( 11 
۰ ۰ 00 ۰2997« 03- 9000۵۰:0۹ 00010 
01 947 ۱ 3 735 20.8090 090106 
00 ۹14 0141 4 1.9881 0352 
00 2-53 ۳2 ۱3 ۳ 4 2.4025 00258 
۵ 00 298« 2 - 5 1.4400 ۱0040 
01 12977 )0 0123- 74 1.5129 0246 
0 ۹214« ۳ 0012 0-4 بووه,و 00530 

00338 4.1209 4 ۳ 003 ۳۰03 00 
)۳23۳07 41209 ۶ - 3 )۳2307 0 . 0 
0490 6.0025 74 - 0.245 0.455 01 
09812 00011 00325 04325 ۳ 
00353 4.4944 7+ «0212 («8212 ۵ 


حال به بررسی طراحی آزمون بهینه می‌پردازيم. هدف از بهینه‌سازی آزمون عمر شتابیده تتش ابت» به دست آوردن مدت‌زمان بهینه اعمال 
آزمون می‌باشد. در این بررسی حامع» عوامل موثر بر زمان بهیته خانمه آزمون از قبیل سطوح تنش» میزان سانسور کردن داده‌ها تعداد 


بازبینی‌ها در هر سطح تنش و حجم نمونه تحت هر دو معیار بهینه‌سازی موردبررسی قرار می‌گیرد. 


موز 
۲ 500۵۲6 ممع]/۱۷ ۶ 
۲ 01211۷۶ 1 3 


جدول ۲- برآوردگر ,۷۲1 پارامترها به همراه اریبی. ۷1517 و 1317 برای داده‌های مثال ۲. 


۰ 62۵۳000۱6 0۲ 0212 ۲0۲ ,۱۸۳ 20 ۱۷۲۵۲ م6 ۱۳۱۸ مصرماه همه تهج 6 ۵۶ ]۱۷۲ -2 016 ] 


۰ 62۵00016 ]0 0202 10۲ ,۱۳ 20 ها م6 ۱۳۱4 مصماه امه هه مظ) ۵۶ ]۱۷ -3 12016 
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"(۳2۰997 
|)۳۰-19956 
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|)۳(2-32(14 
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0۳۰-13 
0141 

039022 
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- (2 
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0121 
0014 
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)(2141 
0۳20022 
-)(.0120 
- 1 
0۳2921 
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1 ۹ 
> 
> 7 
* 2 
1 ۹ 
> 7 
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1 2 
* + 
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(9 
1-3600 
1-401 
1-3600 
210 
625 
1-81 
210 
114100 
7.61 
3022 
1209 


1۳۳ 
)(۰035 
0.00 
0.0600 
9۹99۳9 
0.00520 
0۳125 
0۳2-305۵ 
|) 
۳2-1920 
0.1405 
0.1405 
)۳2-30/7 


جدول ۳- برآوردگر ,]1 پارامترها به همراه اریبی. ۷15 و 1817 برای داده‌های مثال ۳. 


1 1۷۸ و۳12 1۷۹3۳ 1۳1۳ 
0۵ . نم 1.1972 8 .0- ۲ 7.8400 000253 
9 1.8 2 0- 75 4640 .8 08 
80 143 43 .0 75 8490 .1 00039 
0 45 1 45 .0 5 2.0250 0.0028 
0 :0 7 0) 53 ()- ۰:0 :0۵0۳090۹-2 0 (0) 
0۳ 73 .1 7 ,0- 4 6129 1 0009 
0 16 0.6 + 1236 .1 0.0096 
0 52 .1 2 0 75 2.7040 0.032 
۵ 40 1.13 7 0- 4 6129 .1 0.6 
0۳ 1.09 01- 074 941 .5 00176 
0 91.27 0.1 ۷04 33604 0.055 
0 1.23 0.123 +۰ 5129 .1 0.007 


حدول ۹4 مقادیر بهینه زمان خاتمه آزمون را برای داده‌های مثال‌های 3 ۲ و ۳ تحت هر دو معیار بهینه‌سازی و به ازای 6 < 7۲ ارائه 


1 


می‌دهد. در این حدول زمان بهینه خاتمه آزمون تحت معیار اول را با نماد ر77 و زمان بهینه خاتمه آزمون تحت معیار دوم را با نماد 
نمایش می‌دهيم. با توحه به زمان‌های بازرسی از پیش تعیین‌شده برای هر یک از مثال‌هاء نتیجه می‌شود زمان بهینه خاتمه آزمون در تمامی 
مثال‌هاء بلافاصله بعد از آخرین بازرسی قبل از خاتمه آزمون است. بدین ترتیب می‌توان نتیحه گرفت زمان آخرین بازرسی قبل از خاتمه 
آزمون. مزثرترین فاکتور در تعیین زمان بهینه خاتمه آزمون است. به‌علاوه از جدول 5 نتیحه می‌شود که زمان‌های بهینه‌شده تحت هر 
دو معیار یکسان هستند البته همان‌طور که مشاهده می‌نمایید. تابع هدف تحت معیار دوم در هر سه مجموعه داده موردبررسی کمتر از 
مقدار تابع هدف تحت معیار اول می‌باشد» این امر نشان‌دهنده این است که آزمون بهینه طراحی‌شده تحت معیار دوم کاراتر از معیار اول 


4 


است. 





۸۳ 


سم 


براورد بارامترها با استة 


ده از الگوریتم 7و 


بهینه‌کردن طرح | 


زمون عمر شتاییده تنش تاه 


بت با پا 


رامترهای غیر ثا 
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سانسور فزاینده نوع اول 
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تا دوره س شماره ۲ ز هت از ۰-۰۰( صفحه: ۰-۲۱ ۷ ۵ 


جدول ۳- زمان بهینه خاتمه ازمون و تابع هدف متناظر با ان تحت دو معیار بهینه‌سازی. 
۰ ۱]۲۷ هنن و۲ نا مصمتا ص رها فص رومام 20 متا م)مصتصصنع) امصصت)و0 -4 ۵016 ] 


معیار بهینه‌سازی دوم معیار بهینه‌سازی اول 


داهای‌ها.۱ صام تته ات ] من 01 2 7 01 2 :1 

مقارفانع فوف 1 - وه 8 2 رب 
ایشا ۲ اش شاه ابرم 001 2 7 001 2 1 

مقدار تابع هدف 61 -< و 02 2 ره 
تا ان امش تاه مرت 01 2 ۲" 001 2 ,1 

مد تاه قرف 2-8 و۵ 7 2 رره 


زمان بهینه خاتمه آزمون نشان می‌دهد؛ که واضح است. میزان سانسور میانی تأثیری بر زمان بهینه‌شده ندارد. 


جدول ۵- تأثیر احتمال سانسور میانی بر زمان بهینه خاتمه آزمون تحت هر دو معیار بهینه‌سازی. 
۰ 0۵6۵۵۱۲۲ ۲۵ ۲ع0 تا ممحصل) مامتها افصصتامون مط) دم مصتهمعصی ع)ها0عصه ها ۵۴ ۵160 ۰16۵ -5 12016 


7 
0/۵ 10۳2 20/0 20/0 10/0 20/0 600۳/0 0۳/0 50۳/0 90/6 
منال ۱ 1 1 0.3001 0.3001 0.5001 0.5001 0.5001 ۰ 0.5001 ۰ 0.5001 ۰ 0.5001 ۰ 0.5001 
10 1 0.3001 0.5001 0.5001 0.5001 0.5001 ۰ 0.5001 ۰ 0.5001 ۰ 0.5001 ۰ 0.5001 
منال ۲ 1 1 0.300 0.3001 0.3001 0.3001 0.3001 0.3001 0.3001 ۰ 0.3001 ۰ 0.3001 
10 1 0.300 0.3001 0.3001 0.3001 0.3001 0.3001 0.3001 ۰ 0.3001 ۰ 0.3001 
منال ۲ 1 1 0.3001 0.5001 0.5001 0.3001 ۰ 0.5001 ۰ 0.5001 ۰ 0.5001 ۰ 0.5001 ۰ 0.5001 
106 1 0.3001 0.5001 0.5001 0.5001 0.5001 ۰ 0.5001 ۰ 0.5001 ۰ 0.5001 ۰ 0.5001 


شکل‌های ۰۱ ۲ و ۳. نمودار تأثیر احتمال سانسور میانی را بر مقدار تابع هدف بهینه تحت دو معیار موردبررسی به ترتیب برای مجموعه 
مقادیر اولیه مثال‌های ۰۱ ۲ و ۳ نشان می‌دهند. همان‌طور که از نمودارها نتیحه می‌شود. با افزایش مدت‌زمان آزمون» مقدار تابع هدف بهینه 
تحت دو معیار بهینه‌سازی دارای روند افزایشی مشابهی می‌باشد. به‌علاوه مشاهده می‌شود با افزایش احتمال سانسور میانی» توابع هدف در 
هر دو معیار بهینه‌سازی مقادیر کمتری را اختیار می‌کنند» از طرفی هرقدر زمان خاتمه آزمون» کمتر باشد و به سمت صفر میل کند» نیز 
نتیجه می‌شود مقدار تابع هدف کمترین مقدار را به خود می‌گیرد. درواقع می‌توان گفت برای دست‌یابی به آزمون کارا یک راه‌کار بالابردن 
احتمال سانسور میانی واحدها در هر بازرسی است و راه‌کار دیگر کم کردن زمان خاتمه آزمون است. همچنین مشاهده می‌شود که مقدار 
تابع هدف تحت معیار بهینه‌سازی اول به‌طورکلی بیشتر از مقدار متناظر آن تحت معیار بهینه‌سازی دوم است و این بدان معنی است که 
به‌طورکلی آزمون بهینه تحت معیار بهینه‌سازی دوم کاراتر از معیار اول است. 


جدول 1 تأثیر افزایش حجم نمونه را بر زمان بهینه خاتمه آزمون تحت هر دو معیار: به ترتیب برای داده‌های مثال‌های ۰۱ ۲ و ۳ با احتمال 
سانسور میانی ۱۰۵ نشان می‌دهند. در حدول ۲" طی چهار حالت ححم نمونه‌های مختلفی برای مثال‌های ۰۱ ۲ و ۳ در نظر گرفته‌شده 
است. برای مثال ۱: حالت اول 220 ره7< ,7 حالت دوم 0 - م< رل حالت سوم 0 - ره ,7 و حالت چهارم 50 - وه7< ر7. 
برای مثال ۲: حالت اول 220 و7- وه رهه» حالت دوم 30 - > رمرم حالت سوم 40- و- ر< رم و حالت چهارم 
0- و - و2 7 و برای مثال ۳: حالت اول 20-ب- و- و7<ره. حالت دوم 30-"- 7- و7 ,«. حالت سوم 
0 - پا - 91 - ر< 7 و حالت چهارم 0  <‏ - 1 - 7 ,7 . با توحه به داده‌های حدول 1 نتیجه می‌شود که با تغییر حجم نمونه. 


در زمان بهینه خاتمه ازمون تغییری به وحود نمی‌آید. 
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76۲۳۱۱۳۵۵ ۵ 6۲۳۱۱۳۵۵ ۵ 776۲۳۱۱۳۵۱۵۲ ۵ 


شکل ۱- تابع هدف بهینه نسبت به زمان خاتمه آزمون به ازای احتمال سانسورهای میانی مختلف. برای داده‌های مثال ۰۱ 
۲ 00050۲1۳۵ ۱۱0۵)۵۱۵۵۵60126 )00۵ ۲۵ ممصصرن) 0)هصتصصه ما هرمع ۱۲۱6 0ص ماع زهان اهصصتان0۵0۵ -1 عاعز۲ 
۰ 62100016 0 0212 


سم 


2500 


م72 


:2 0 


بهینه‌کردن طرح | 


سح 





۲۳۱۱۳۵۵ ۵ 7۲6۲۳۱۱۳۵0۵ ۵ 26۲۳۱۱۳۲۵۵ ۵ 


شکل ۲- تابع هدف بهینه نست به زمان خانمه آزمون به ازای احتمال سانسورهای میانی مختلف. برای داده‌های مثال ۲. 
۲ 00050۲1۳0۵ ۱۱۵6۵۱060126 )00۵6 ۲۵ ممصصرن) ص0)هصتصصه م6 هرمع ۱۲۱6 0ص ۲ تععزهان افصصتا0۵۵ -2 ۲1۵۱۵۲۵ 
۰ 62100016 0 0212 


زمون عمر شتابیده تتش ثا 


بت 


7۳-26 


با پارامترهای غیر ثا 
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4 مه 


0 05 ۱ 05 
7۲6۲۲۱۱۳۵۵ ۵ 76۲۲۱۱۳۵۵ ۵ 7۲6۲۳۳۱۳۵۹۵ ۵ 





سانسور فزاینده نوع اول 


شکل ۳- تابع هدف بهینه نست به زمان خانمه آزمون به ازای احتمال سانسورهای میانی مختلف. برای داده‌های مثال ۲. 
۲۲ 000501108 ۵۱0۵601266) ۱۱ )ص006 ۲۵ ممصصر) 0)هصصع ما هرمع ۲۱۱6 ۵ص ماع زهان افحصتام۵ن0 -3 ۲۵0۲۵ 
۰ 200016 0 0219 


جدول 1 - تأثیر حجم نمونه بر زمان بهینه خاتمه آزمون تحت هر دو معیار بهینه‌سازی. 
۰ هنن ۲۵ مها مممصا) م60 مضه مط چم م62 عاصرصصنهی ۵ )66 1۳6" -6 12۱016 





مثال ۱ 0 - ۸۷ 0 - ۷ 6۱0 - ۷ 0 - ۷ 
,1 0(2-301) 02-001) 0(2-001) 0001) 
106 0(2-301) 0(2-2301) 0(2-01) 05001) 
مثال ۲ 0 - ۸ 0 - ۸ 0 - ۸ 0 - ۸ 
,1 03001) 3001() 0.3001) 001(-0(2-3) 
16 03001) 03001) 0.3001) 0(2-32-001) 
منال ۳ 0 - ۸۷ 0 - ۸ 120 - ۷ ۵0 - ۷ 
,1 0(2-01) 0(2-2301) 0(2001) 0001) 
16 0(2-232001) 0(2-2301) 0(2-3001) 0001) 
9۸۷۹ 
در شکل‌های 6 تا ۰٩‏ نمودارهای تأثیر افزايش حجم نمونه بر تابع هدف بهینه متناظر با جدول " برای سه مثال موردبررسی ترسیم‌شده 
دست‌یابی به آزمون بهینه کارا بهتر است از حجم نمونه کمتر استفاده شود. 
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شکل ۴- تأثیر حجم نمونه بر تابع هدف بهینه تحت هر دو معیار برای داده‌های مثال 1 


۰ 62100016 0۴ 0۵1۵ ۲۵۲ راهان )۵ م0 0ص ۲۵ رصان اهصصتامرن عظ) مه 6۱2۵ عامرصصوی 01 601160 116 -4 ۲۱۵۱۵۲۵ 
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400 ِ 


یم‌گیری و تحقیق در عملیات. دوره ۶ شماره ۰۴ زمستان ۱۴۰۰ صفحه: ۵۷۰-۵۹۱ 
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شکل ۵- تأثیر حجم نمونه بر بر تابع هدف بهینه تحت هر دو معیار برای داده‌های مثال ۲. 
۰ 6۵100۱6 0۴ 0212 0۲ راهان ۵۵ 0ص ۲۵ رها افحصتامون معط ده م62 عامرصصوی 0۴ 61160 11۳6 -5 ۲۱۵۱۲6 
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شکل 1" - تأثیر حجم نمونه بر تابع هدف بهینه تحت هر دو معیار: برای داده‌های مثال ۳. 
۰ ۵۵10016 0۴ 0242 ۲۵۲ ردتهاتن )۵۵ م0صنا متا زان افصهتامن مط من 5126 مامرصصهه ]0 ۵11606 1۳۵ -6 ۲۱۵۱۵۲۵ 


حدول ۰۷ تأثیر افزايش تعداد سطوح تش را بر زمان بهینه خاتمه آزمون تحت هر دو معیار موردبررسی برای داده‌های مثال‌های ۰۱ ۲ و ۲ نمایش 
می‌دهد. در تمامی این مثال‌ها؛ تعداد بازرسی 1 و زمان‌های بازرسی 9 2 ۶« هِِ« و ۵ در نظر گرفته‌شده ادن همان طور 
که مشاهده می‌شود. زمان بهبه خاتمه آزمون متأثر از تعداد سطوح تنش نیست. 


جدول ۷- تاثیر تعداد سطوح تنش بر زمان بهینه خاتمه آزمون تحت هر دو معیار بهینه‌سازی. 
۰ ۱۲۷ هنن 6۲۵ عص ممصص ملامصصع) ماع مرن مزم۲ع۱ عممای ۵ اجه عظع ۵۴ عم 16 -7 12016 


4 - 1 ۲-3 2 - ۲ 
مثال ۱ . م1 05001 0500 0۹00 ۰ 
1 05001 05001 05001 
مثال ۲ ِ 05001 05001 05001 
ِ 05001 05001 05001 
منال ۳ ۳ 05001 05001 05001 
1 05001 05001 05001 


شکل ۷ نمودار تأثیر افزایش تعداد سطوح تنش را بر تابع هدف بهینه با توحه به هر دو معیار موردبررسی برای داده‌های مثال ۰۱ نمایش 
می‌دهد. همان‌طور که مشاهده می‌شود با افزایش تعداد سطوح تنش» تابع هدف در هر دو معیار افزایش می‌یابد. پس برای دست‌یابی به 


آزمون بهینه کارا بهتر است از تعداد سطوح تنش کمتری استفاده نماییم. 
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شکل ۷- تأثیر تعداد سطوح تنش بر تابع هدف بهینه تحت هر دو معیار. 
۱ ۱6۲ص ۲ رصان افصانن معط من ولع۲ع۱ ماو 0 ماه ان 61160 116 -7 ۲۱۵۱۵۲۰۵ 
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سانسور فزاینده نوع اول 


حدول ۸ تأثیر افزایش تعداد بازبینی را بر زمان بهینه خاتمه آزمون تحت هر دو معیار موردبررسی برای داده‌های مثال‌های ۰۱ ۲ و ۳ نشان 
می‌دهد. در این حدول برای 1-2 زمان‌های بازرسی را ۰۰,۱ ۸/۲ برای 1-3 زمان‌های بازرسی را ۰۰,۱ ۰۰,۲ ۰,۳ برای 1۱-4 
زمان‌های بازرسی را ۰۰۱ ۰۰,۲ ۰۰/۳ ۰4 برای 1<3زمان‌های بازرسی را ۱ ۰۰ ۰۰,۲ ۰۰,۳ ۰۰,4 ۰/۵ وبرای 1-6 زمان‌های 
بازرسی را ۰۰7۱ ۰,۲ ۰۰/۰۰/۳ ۰/۸۰۵ در نظر گرفته‌شده است. با توجه به نتایج جدول ۸ متوجه می‌شویم که زمان بهینه خاتمه 
آزمون تحت هر دو معیار بلافاصله بعد از زمان بازرسی قبل از خاتمه آزمون است. 





جدول ۸- تاثبر تعداد بازبینی بر زمان بهینه خاتمه آزمون تحت هر دو معیار بهینه‌سازی. 
۰ هنن 6۲۵0 ۱06 ممصصت لته مظن مصرمتایءعن مه ۵۶ مج صصناهد ماع ۵ 01160 116 -8 2016 ] 


۸۸ 6 - 111 5 - 17 4 - 1 3 - 1۳ 2 -7 17 
مثال ۱ ,1 1() 0-1) 1) 0001) 01001) 

)01001 )0001 )0(01 )0-1 )(1 10 

منال ۲ ,1 1() 0-1) ۵(1) 0001) 01001) 

)01001 )0001 )01 )0-1 )(1 1 

فتالن ۲ ,1 1() 0-1) 0(1) 01) 01001) 

)01001 )001 )001 )0-1 )(1 10 .< 


شکل . نمودار تأثیر افزایش تعداد بازبینی را بر تابع هدف بهینه تحت هر دو معیار بهینه‌سازی برای داده‌های مثال ۱ نشان می‌دهد. ازآنحاکه 


حکمی‌پور / ت 


در بررسی مثال‌های 3 ۲ و ۲ همواره نتایج مشابهی حاصل‌شده و طبق بررسی‌های انجام‌شده تعداد سطوح تنش تأثیری درروند تعبین زمان 
بهینه خاتمه آزمون و تابع هدف بهینه ندارد و از طرف دیگر به دلیل محدودیت صفحات مقاله» برای این بررسی فقط از داده‌های مثال ۱ 





۷ ۳ همان‌طور که مشاهده می‌ شود افزایش تعداد ثار ی منحر به افزایش میزان تابع هدف تحت هر دو معیار بهینه‌سازی 
۳ 5 ِ 
۳ می‌شود. درنتیحه برای دست‌یابی به آزمون بهینه کارا بهتر است تعداد بازبینی‌ها در حین آزمون کم انتخاب شود. همچنین مشابه نتایحی که 
۳ از شکل ۱ تااشکل ۷ حاصل شد. در شکل ۸ نیز واضح است که معبار بهینه‌سازی دوم نسبت به معیار اول» کاراثر است. 
2 
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شکل ۸- تأثیر تعداد بازرسی‌ها بر تابع هدف بهینه تحت هر دو معیار برای داده‌های مثال ۱. 
8 ۲0۲ ,هن ۵۵6 عصنا ماع زهان امصتامون مط دص ممصص انم 0۶ رم صاحصناه ۵ 61160 11۳6 -6 ۲۱۵۱۵۲۵ 
1۰ 200016 0۴ 


۶ مثال کاربردی 


و این بخش.» موضوعات بررسی‌شده در قسمت‌های قبل را برای یک محموعه داده واقعی. به کار می‌بریم. داده‌ها برگرفته از کتاب تلسون 
( است. در فصل سوم این کتاب نتایج حاصل از آزمون عمر شتابیده تتش ثابت برای دستیابی به زمان خرابی عایق ترانسفورماتور 
ارائه‌شده است. این آزمایش شامل اعمال ولتاژ در سه سطح ثابت و بیش ‌ازحد طبیعی می‌باشد که به ترتیب عبارت‌اند از ۳۵,۴ کیلووات 
۴ کیلووات و ۴۶/۷ کیلووات. در هر سطح تنش ۱۰ واحد آزمونی تحت آزمایش قرارمی‌گيرند. لازم به ذکر است» سطح ولتاژ در 
حالت استفاده معمول در این آزمون ۱۴,۴ کیلووات است. داده‌های این مثال در حدول ٩‏ ارائه‌شده است. برای این داده‌ها توزیع لومکس 


را با پارامترهای *2.1683<10 ۵ و ۸24.909 برازش می‌دهيم. با توحه به اينکه در آزمون کلموگروف-اسمیرنوف" مقدار آماره 
2021 1-۵ و 0.5/43< 7-۷6 است. می‌توان نتیجه گرفت این برازش, برای داده‌ها در شرایط استفاده معمول» مناسب است. 
اکنون ما بر روی‌داده‌های این مثال» روش بازرسی دوره‌ای و سانسور فزاینده نوع اول را اعمال می‌کنیم. برای بازرسی دوره‌ای» دو زمان از 
پیش تعبین‌شده ۲۰ و ۴۰ را به‌عنوان زمان‌های بازرسی در نظر می‌گیریم که سانسور فزاینده در همین زمان‌ها با احتمال سانسورهای میانی 
معادل 6۰ 50۱۰ و 50۲۰ به واحدهای تحت آزمون اعمال می‌شود. با استفاده از الگوریتم 7711 برآورد ماکسیمم درستنمایی پارامترهای 
محهول مدل» تحت هر یک از سانسورهای میانی ذکرشده به همراه خطای استاندارد" (92)متناظر آن‌ها در حدول ۱۰ ارائه‌شده است. 
در حدول س ۱ مشاهده می‌شود که خطای استاندارد حاصل از بر آورد با استفاده از الگوریتم 171 مقدار کمی است که این تسان از 
مناسب بودن روش برآورد بررسی‌شده در این مقاله دارد. به‌علاوه مشاهده می‌شود که با افزایش سانسور میانی» مقدار خطای استاندارد به 


میزان کمی افزایش می‌یابد و این با نتبحه حاصل از مطالعات شبیه‌سازی که در بخش قبل ارائه شد منطبق است. 


جدول -٩‏ زمان خرابی عایق ترانسفغورماتور. 
۰ 0۳20601۱۵۱6۱۰۶ ۵۴ ۱۵۵5 ۲۵۱[۵۲6 -9 عوج1 


ولتاژ زمان‌های خرابی 

۴ کیلووات 01 4 5 .۳ 2047 7 ۰ 3083 53739 10023 1002.3 
۴ لو وان 06 134 2 ۵ 199 250 2 . 328 4444 502 56.2 
۷ کیلووات 31 83 و 90 136 149 ۰-61 169 213 481 


جدول ۱۰ - برآوردگر ,]۷۲ پارامترها به همراه 917 برای داده‌های واقعی. 
۰ ۲۵۵۱ ۲0۲ ره ارعط) ۱۱۱۵ مجماه همامه هن عظ) ]۵ بزر[۱۷ -10 12016 


ِ 00 1" 00 01 
0۳/6 ۳۲ ]1۷ 00043) 107-.۳(2) 1-04 2.4 
۳ 00295) 0.1959 3069 -0(2) 00193 
۵ ار( ۳2."(0(037) 02073 1-1206 2.045 
۳ 03753) 0۳1977 9069(-۳(2) 00225 
200۳6 00 9039(-۳(2) 0.1973 1-3992 19975 
۳ ۳23929) 0.2195 0200935) 0.235 


در حدول ۱۱ مقدار بهینه زمان خاتمه آزمون تحت سانسور میانی ۲۰۹۵ ارانه‌شده است. همان‌طور که مشاهده می‌نمایید» مشابه نتایحی 
که از مطالعات شبیه‌سازی گرفته شد. زمان بهینه خاتمه آزمون در هر دو معیار باهم برابر و بلافاصله بعد از زمان بازرسی اول می‌باشد. 


همچنین مقدار تابع هدف متناظر با معیار بهینه‌سازی دوم کمتر از مقدار تابع هدف متناظر با معیار بهینه‌سازی اول است. 


جدول ۱۱- زمان بهینه خاتمه آزمون و تابع هدف متناظر با آن تحت دو معیار بهینه‌سازی برای داده‌های واقعی. 
۲ ه0۵1 ۲م اون و ما مصمتاص زان مصی واه 220 ممصص 10)هتصصع) افصصتان0۵۵ -11 12016 








۰ 1۲62۸۱ 
معیار بهینه‌سازی دوم معیار بهینه‌سازی اول 
زمان بهینه خاتمه آزمون 1 -< 6[ 1 0 << ]1 
مقدار تابع هدف 40 < و 7 -< ,6 


۷- نتیجه‌گیری 
در این مقاله آزمون عمر شتابیده تنش ثابت ] مرحله‌ای را برای محصولاتی که طول عمر آن‌ها از توزیع لومکین. غیت هی کین به 


کاربردیم. همچنین برای کاهش هزینه‌های آزمون از سانسور فراینده نوع اول و روش بازبینی دوره‌ای برای جمع‌آوری داده‌ها استفاده 


۷مصنصک 010۷عمصصام 1 
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یم‌گیری و تحقیق در 


مب 


تا دوره س شماره ۲ هت از ۰-۰۰«( صفحه: ۰-۱ ۷ ۵ 


کردیم. برای تطابق بیشتر آزمون با آنچه در عمل اتفاق میفتد» در نظر گرفتیم که تتش اعمال‌شده بر هر دو پارامتر توزیع لومکس تأثبرگذار 
است و برای بررسی تأثیر افزایش سطح تدش بر طول عمر محصولات. فرض شد که هر دو پارامتر مقیاس و شکل دارای رابطه لگ خطی 
هستند. در این مقاله برای برآورد ماکسیمم درستنمایی پارامترهای مجهول از الگوریتم بهره گرفتيم. اصل اطلاعات ازدست‌رفته به 
الگوریتم 7711 کمک می‌کند تا اطلاعات موحود در داده‌های مفقودشده را به کار گیرد. عملکرد این روش برآورد از طریق معیارهای اریبی» 
و 1/5 مورد ارزیابی قرار گرفت. در ادامه مقاله به‌منظور بهینه‌سازی طرح این آزمون از دو معیار بهینه‌سازی مبتنی بر ماتریس اطلاع 
فیشر استفاده کردیم و مدت‌زمان بهینه اعمال هر سطح تدش را با توحه به هر دو معیار محاسبه کردیم. سپس تأثیر فاکتورهایی از قبیل حجم 
نمونه. نسبت سانسور میانی» تعداد سطوح تنش و تعداد بازبینی بر مدت‌زمان بهینه‌شده و کارایی نسبی آزمون توسط مطالعات شبیه‌سازی 
شده و یک مجموعه داده واقعی موردبررسی قرار گرفت. با توجه به بررسی‌های صورت گرفته مشاهده می‌شود. زمان بهینه‌شده در تمامی 
مثال‌ها تحت هر دو معیار یکسان است و این زمان بهینه بلافاصله بعد از آخرین بازرسی قبل از خاتمه آزمون می‌باشد. بدین ترتیب می‌توان 
نتیجه گرفت زمان آخرین بازرسی قبل از خاتمه آزمون» مزثرترین فاکتور در تعبین زمان بهینه خاتمه آزمون است. به‌علاوه تحت معیار 
بهینه‌سازی دوم. آزمون کاراتری نسبت به معیار اول حاصل خواهد شد. همچنین برای دست‌یابی به آزمون کاراتر بهتر است از حجم نمونه 
کم تعداد سطوح تنکن کمتزه تعداد باز نیش کمتر و اختمال ساتسور میات بنتتر استفاده کرد به‌این ترتیتب:هی توان ادعا کرد که‌به آزمون 
کارا دست‌یافته‌ايم. 
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9۹۱ 


سم 


براورد پارامترها با استة 


ده از الگوریتم 7و 


سح 


بهینه‌کردن طرح آزمون عمر شتابیده تتش ثا 


بت 


« 


با پا 


رامترهای غیر ثا: 


بت تحت 
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